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Abstract 

We present an experimental investigation of a novel low Reynolds number shear flow instability 
triggered by a chemical reaction. An acid-base reaction taking place at the interface between a 
Newtonian fluid and Carbopol-940 solution leads to a strong viscosity stratification, which locally 
destabilizes the flow. Our experimental observations are made in the context of a miscible dis- 
placement flow, for which the flow instability promotes local mixing and subsequently improves 
the displacement efficiency. The experimental study is complemented by a simplified normal mode 
analysis to shed light on the origin of the instability. 
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I. INTRODUCTION 



Designing a method to locally control the hydrodynamic stability of a shear flow is of 
importance in many practical and laboratory applications. In the absence of inertia, non- 
equilibrium hydrodynamical systems may lose stability if a relevant physical field becomes 
stratified. To help illustrate this point, convective motion may be induced in a thin static 
layer of fluid when heated from below [l| or gravity induced density stratification may sustain 
internal gravity waves, 2|. An alternative method of locally destabilizing low Reynolds 
number (Re) shear flow would be to induce significant changes in the local fluid rheology. 
Although flows with a strongly stratified viscosity have been proved to be theoretically 
unstable, this situation has never been demonstrated experimentally. With flows of 
simple Newtonian fluids, it is difficult to vary the viscosity locally to induce an instability. 
With complex (or structured) fluids , however, the situation is significantly different: the 
rheology is strongly coupled to the molecular scale organization of the fluid. This opens a 
new possibility of locally controlling the viscosity by inducing local changes in the molecular 
structure via a chemical reaction. The advantage of such method is that a chemical reaction 
may be controlled by either mass transfer or by local heating or cooling. In this paper we 
show experimentally that a displacement flow of two miscible liquids may be destabilized 
by local changes in the fluid rheology triggered by an acid-base reaction at their interface. 

Miscible displacements have been studied in depth within the context of Hele-Shaw and 
porous media displacement instabilities, initially motivated by a desire to better understand 
oil reservoir recovery issues, e.g. 

aaaaaa 

, and this science is now well developed. 
In the Navier-Stokes setting, miscible displacements through small ducts have been studied 
both experimentally and numerically, 10, ll|, Q|- For the large Peclet number (Pe) 
regime, depending on the viscosity ratio m, quasi-steady viscous fingers may form and prop- 
agate with a sharp displacement front that is retained over long timescales. The efficiency 
of the displacement (defined as the amount of fluid removed from the pipe walls) may be 
either measured or estimated using computation or asymptotic methods, see e.g. [3, Q. 
In the case that the fingers remain stable, the residual wall layers thin at a rate that is sig- 
nificantly slower than the mean flow and the efficiency of the displacement is < 1. Although 
asymptotically the efficiency may approach 1 as t — > oo, for various practical reasons one 
might not want to wait. 
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To improve the displacement efficiency, one method would be to trigger a flow instability, 
such that the perturbed interface disturbs the residual wall layers. At large Re, multi-layer 
flows of viscous fluids are usually unstable, but linear interfacial instabilities are also found 
for quite low Re. Dating from the late 1960's, there are a number of studies involving both 
immiscible and miscible fluids, e.g. {15I . An extensive review can be found in the text 



171 ]. and the physical mechanisms governing short and long wavelength instabilities have 
been explained by [3, [if]]. These studies generally refer to the situation where there is a 
jump in the viscosity at the interface between two fluids. If the change in viscosity is instead 
gradual, e.g. due to a diffuse interfacial region, then the stability characteristics appear to 
mimic those of the system with the viscosity jump, see 3|. 

In porous media, non-monotone viscosity variations are known to cause linear instabil- 
ity of planar displacements, see j^, Q, 21], as can be predicted by classical mobility ratio 



arguments. Here however, we are in the Navier-Stokes regime and consider primarily shear 
flows. Also we have sharp localized change in viscosity, which is hard to achieve in "sim- 
ple" fluids, where viscosity is often related to slowly varying concentration or temperature, 
i.e. due to molecular structure of the fluid. For complex (or structured) fluids the vis- 
cosity depends strongly on the microscopic structure. In the case of polymer solutions, 
the microscopic structure of the fluid can be locally modified by either mechanical means 



(e.g. shear-thinning, 



22J) or chemical means, by locally modifying the chemical bonds be- 



tween neighboring polymer molecules. 

In recent years there have been a number of studies of systems with coupled chemi- 
cal reactions and fluid flow, so it is natural to examine the relation to this literature. In 
the first place even without fluid flow, spatially traveling waves can be observed in chem- 
ically excitable media governed by coupled reaction-diffusion systems, see e.g. the review 
23| . These fronts may frequently destabilize linearly. There is often coupling between the 
reactions, and hence a mechanism of feedback, and sometimes a significant difference in 
the reaction rate constants. Slightly simpler systems involve single species auto-catalytic 
reactions, which mathematically admit traveling (chemical) wave solutions. Systems in 
which the auto-catalytic reaction results in a significant density change have been studied 
intensively. The base traveling chemical wave is coupled with a Rayleigh- Taylor problem. 
There are consequently a wide range of stable and unstable situations to be studied, and 
even more once a base fluid velocity is considered. A selection of the many works includes 
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3JJ. Of closer relation to our work is the sequence of papers, 



351 ] . which have considered a concentration dependent viscosity in the context of 
miscible porous media displacements. Although interesting and a useful guide in method- 
ology, direct relevance of these many papers to our study is in fact limited. Our system is 
not auto-catalytic, there are no buoyancy effects and the displacement flow is not a gradient 
flow. 



A. Industrial motivation 

The motivation for our study comes from the construction of oil and gas wells. Since the 
early 1990's there has been an increasing number of wells that are constructed with long 
horizontal sections. The worlds longest extended reach wells have horizontal sections in the 
10 — 15 km range, but these are exceptional. More routinely, wells are built with horizontal 
extensions of up to ~ 7 km. One of the key barriers in constructing longer wells comes from 
simple hydraulic friction. In a vertical well, both the pore pressure of reservoir fluids and the 
fracture pressure of the reservoir rock increase with depth, approximately linearly. Judicious 
choice of fluid density and circulating flow rates keeps the wellbore pressure inside the so- 
called "pore-frac envelope", i.e. the region where the porous rock does not fracture. In a 
horizontal well section, the pore-frac envelope is unchanged with length along the well, but 



the 
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rictional pressure increases with length, leading to eventual breaching of the envelope 
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38]. 



There is a consequent interest in methods and fluids that control the frictional pressure 
in some way. Two operations where this is important are drilling and cementing. In the 



cementing process, [39j], it is necessary to displace the drilling mud with a spacer fluid and 
then with a cement slurry. As the section is horizontal, density differences between the fluids 
lead to stratification and should be avoided. Instead the focus is on controlling the rheology 
of the fluids and the displacement flow itself. The idea behind the reactive instability that 
we are studying is explained by the following simple calculation. 

Suppose simplistically that we have Newtonian fluids, a circular pipe of radius R and that 
we wish to displace at mean speed Uq. If Ap e is the difference between pore and fracture 
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limits, then we are restricted to a length-viscosity combination: 

R 2 

Normally, we will need to select the viscosity of the spacer fluid, fj s > mfjd, where m > 1 
and fjd is the viscosity of the in-situ fluid, e.g. drilling mud, so that 

(1) 

8mfj d U 

Suppose instead we are able to pump a spacer fluid of viscosity fj s < fj d , that reacts with the 
in-situ fluid to cause a local instability. If the instability causes effective mixing across the 
pipe after the fluid has traveled k radii, we may model this process by a radial diffusivity 
D r ~ (RU )/k. Provided k <C L/R, we have an effective Taylor dispersion process in 
which the mixed zone diffuses axially along the pipe, relative to the mean flow. This axial 
diffusion of the mixed region is governed by a dispersion coefficient of order ~ kDT{RUo), 
where Dt = 1/48 is the Taylor dispersion coefficient for a pipe. After traveling a distance 
L, the mixed zone will have dispersed axially a distance 

1/2 



L m i x ~ 



kD T RL 



If the reacted mixture has viscosity (a+ l)rj s , where a > 0, then the frictional pressure drop 
along the pipe can be limited by 

, y s Uo 
R 2 

and the pipe length restriction is 



Ap e > 8^ (L + aL r 



L < , R ^ , , . (2) 
^U^l+alkDrR/L} 1 / 2 ) 

Therefore, for a given pore-frac limit, we may increase the length of the well that may be 
effectively displaced provided that: 

^ ; . > 1. (3) 

f) s {\ + a[kD T R/ 'L] 1 /*) 

Over lengths sufficiently long that a[kDTR/ L] 1 ^ 2 1, we achieve a length increase by a 
modest factor mfjd/fjs- 

Even modest increases in length may have a large impact, both economically and envi- 
ronmentally. In modern offshore drilling, one vertical well drilled down from the seabed may 



act as the stem for many lateral horizontal branches, extending radially outwards. Thus, 
increases in the length of horizontal branches correspond to (length) 2 increases in the area 
of reservoir that may be reached and contribute to a reduction in the number of wellheads 
required. This reduces the cost of the field development reduces the environmental foot- 
print, lessens risks of leakage of reservoir fluids and eventually makes for an easier well 
abandonment. 

B. Organization 

Our paper is organized as follows. Section [IT] describes the experimental setup, the mea- 
surement techniques and the rheology of our fluids. Experimental results are presented in 
§1111 Section IIVI introduces a simple hydrodynamic stability model that gives insight into 
the instabilities we observe. The paper closes with a brief discussion of our findings and a 
discussion on future theoretical and experimental studies. 

II. DESCRIPTION OF THE EXPERIMENTS 
A. Experimental apparatus and techniques 

All experiments were conducted in the apparatus illustrated schematically in Fig. HJ It 
consists of a horizontal flow channel FC with circular cross section of radius R = 5.5 mm and 
length L = 1.2 m, immersed in a water filled glass container GC to ensure distortion-free 
flow illumination and imaging. The system was illuminated by a thin laser sheet passing 
horizontally through the transparent walls of both the glass container and the flow channel, 
at the middle vertical position. The laser sheet has a thickness of approximately 50 /im in 
the center of the set-up and about 160 /^m near the walls of the channel. It was generated 
by passing a laser beam delivered by a 150 mW solid state laser, L, through a block of two 
crossed cylindrical lenses, CO, mounted in a telescopic arrangement. 

The flow was imaged from the top (Fig. [TJ) with a charge coupled device camera, CCD, 
equipped with a 35 mm photographic lens. The images were digitized with 8 bit quantization 
and 640 x 760 pixels resolution (which accounts for 130 /im). The size of the imaged area 
was thus 1.4 x 1.4 cm 2 . 
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The flow in the channel was induced by a syringe pump SP actuated by a precise stepping 
motor, Mi (from CBVL, Vancouver) and controlled by computer via a serial port. The 
inflow mean velocity was controlled with an accuracy better than 100 fim/s. In order to 
monitor the evolution of the fluid interface downstream, the CCD and the cylindrical optics 
block, CO, were mounted on a linear translational stage parallel to the flow channel (TS 
in Fig. [1]). The stage was actuated by a stepping motor Mi and controlled by computer 
via a serial port. The typical configuration of fluids during the horizontal displacement 
experiments is schematically illustrated in the inset of Fig. [TJ 

The stability of the interface between fluids was investigated using either Laser Induced 
Fluorescence (LIF) , or Digital Particle Image Velocimetry (DP IV) . The image acquisition 
software was developed in-house and allowed us to adjust the time delay between successive 
frames in relation to the local flow velocity (in order to keep the mean particle displacement 
in the range 5 — 30 pixels). For low values of the flow velocity the time delay was 66 ms. For 
higher flow speeds, the delay was decreased to 16 ms. The fluids were seeded with either a 
small amount (approx. 160 ppm) of 20 /im Polyamide spheres (from Dantec Inc.) for DPIV 
measurements, or with fluorescein sodium salt (Sigma Aldrich) for the LIF measurements. 



Time series of the velocity fields were obtained by a multi-pass DPIV algorithm [40J. The 
spatial resolution was 167 /im. The accuracy of the method has been carefully checked by 
running test measurements for low Reynolds number Poiseuille flows with Newtonian fluids 
and comparing then with the analytical solution. 

In addition to these techniques, we have developed a third visualization method to mea- 
sure the local pH in the flow field. We do so by using a pH sensitive colored dye, Bromoth- 
ymol Blue (Sigma Aldrich). The local value of the pH near the interface was assessed by 
measurements of the color distribution in the field of view. With this we are able to estimate 
the local rheological properties. 



B. Fluid properties 

We have used the same base fluids for all our experiments, with various adjustments to the 
levels of acid and base used in each fluid to examine different regimes. The displacing fluid, 
Fluid 1, was a 65% aqueous sucrose solution. The pH of this fluid was varied between 7 and 11 
by titration with different amounts of NaOH ranging from 150 parts per million (ppm) and 
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300 ppm. The displaced fluid, Fluid 2, was a mixture of a 5% (wt) Carbopol 940 (C-940) 
solution and a 65% (wt) sucrose solution in deionized water. C-940 is generally referred to 
as a weak polyacrylic acid and dissociates in solution. The pH of this solution is typically 
between two and three. During our experiments, neutralization of C-940 molecules present 
in Fluid 2 occurs locally, in the vicinity of the interface, thro ugh contact with Fluid 1. 

The pH of the fluids were measured with a thermally corrected digital pH meter (Topac 
Instruments) with 1.5% accuracy. The rheological properties of both fluids 1 and 2 were 
measured with a stress controlled rotational rheometer: CVOR acquired from Bohlin In- 
struments (Malvern Inc., www.malvern.co.uk). Example flow curves at 22°C of Fluid 1, 
Fluid 2 (at pH « 3) and Fluid 2 (at pH rs 7), are presented in Fig. [2J As shown, the 
viscosity of Fluid 1 is independent of shear rate. With Fluid 2 at a pH m 3.2, the viscosity 
is also approximately constant. As a guide, at a shear rate of 2 s _1 we have viscosities 
fji = 114 mPas and 172 = 138 mPas, for fluids 1 and 2 respectively. As the pH is increased 
the rheology of Fluid 2 changes dramatically (see Fig. [2]). A detailed study of the coupling 
between the pH, the molecular structure and the rheological properties of C-940 has been 



recently reported in 



4l|. 



In order to quantify the dependence of the rheological properties of Fluid 2 on the pH , 
several batches have been prepared at different pH values, and the rheological tests have been 
conducted. The pH dependence of the shear viscosity of fluid 2 at rate of strain 7 = 2 
is shown in Fig. [HK In the range of pH shown, the shear viscosity increases monotonically, 
up to a value nearly two orders of magnitude larger than the non-neutralized value. The pH 
dependence of the measured yield stress of Fluid 2 is shown in Fig. [3b- A further increase of 
the pH (data not shown in Fig. [3]) results in destruction of the gel structure, which causes 
both viscosity and yield stress to drop to their initial low pH values. 



C. Experimental procedure 

Our displacement experiments were conducted as follows. First, the alignment of the laser 
sheet and focus of the camera is carefully checked at several locations downstream. Second, 
the valves V\, V2 are open and the channel is initially filled with Fluid 2. Next, we close 
valve Vx, open valve V3 and start to drive the syringe pump at very low speeds for several 
tens of seconds. This procedure allows us to eliminate any air bubbles generated during the 
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filling of the flow channel. Finally, we close valve Vi and open valves V\ , and set the fluids 
in motion by operating the syringe pump at the desired speed. The flow images were usually 
acquired halfway down length of the pipe, though in some of the experiments the camera and 
the laser sheet were moved at constant speed along the flow channel. Approximately 1000 
flow images were acquired both before the entrance of the fluid interface in the measuring 
window and after. Details of the experimental conditions tested are given in Table [H 

III. EXPERIMENTAL RESULTS 

We propose to control the stability of an interface between a Newtonian and a C-940 
solution by inducing changes in the local fluid rheology via an acid-base reaction, which 
depends on the free charge mismatch across the interface, and the pH of each fluid. With 
1% NaOH, the saccharose solution (Fluid 1) has a pH around 11. When in contact with 
Fluid 2 at pH w 3, neutralization occurs near the interface resulting in a thin layer of high 
viscosity fluid, see Fig. [2j The high viscosity interfacial layer apparently undergoes a self- 
sustained hydrodynamic instability that results in local mixing of the two fluids. Below we 
will present a full description of this reactive instability and displacement. For comparison, 
we have conducted 2 experimental test sequences (at varying flow rates) in which there is 
no reaction, see Table [I] series 1 and 2. In addition we have performed a number of reactive 
experiments, see series 3 in Table [H In total, we performed 5 experiments in this series. 

A. Control experiments 

The first control sequence concerns the two base fluids without reaction. We have seen 
in Fig. [2] that the viscosity of Fluid 2 at pH w 3 is largely constant and governed by 
that of the underlying saccharose solution. We therefore conduct a control displacement 
of Fluid 1, a 65% saccharose solution, by Fluid 2, a 66% saccharose solution. The flow 
rates are chosen so that Re < 1 in all cases, and often Re <C 1. Both fluids are at neutral 
pH. The displacing fluid is mildly more viscous, but not sufficiently so to ensure a good 
displacement, (fji = 114 mPas and r) 2 = 138 mPas) as illustrated in Figure HI The interface 
develops into a long finger that stretches progressively along the pipe. At no time do 
we observe any interfacial instability, and the interface remains sharp, i.e. the length of the 
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tube is insufficient for molecular diffusion to have any significant effect. Qualitatively similar 
results for miscible Newtonian displacements can be found in [lo| . 

The second control sequence also concerns displacement of the two fluids at neutral 
pH, i.e. without chemical reaction, but Fluid 2 has a much larger viscosity (7)2 ~ 20 Pas at 
7 = 1 s~ x ) and a significant yield stress (a y ~ 6.6 Pa). The viscosity of Fluid 1 is unchanged. 



Similar flows have been studied experimentally before in slightly larger tubes 



42 



commented in 



43|. As 



42| . some care is needed in choosing the displacing flow rate, since at low 



flow rates the displacing fluid effectively fractures through the elastic gel, leaving a rough 
edge. For low flow rates we also observed a slightly granular interfacial texture and also 
some asymmetry of the finger. The slight asymmetry of the finger visible in Fig. could be 
due to the non-symmetric entrance condition (the interface between fluids enters the flow 
channel via a T shaped junction, as shown in Fig. [Q which is preserved at all later times 
due to the significant yield stress of Fluid 2. Alternatively the asymmetry may relate to 
a transition between solid-like and fluid-like behavior. We do not attempt to explain this 
further. 

Sample images of an experiment in this sequence are shown in Fig. [5j The C-940 solution 
yields in the center of the pipe, where there is a two-dimensional flow, but the stresses are 
not sufficient to make it yield at the wall. Consequently a static residual wall layer is left in 
the tube as a finger of Fluid 1 penetrates steadily along the pipe. The shape of the nose of 
the finger is quite rounded and the wall layers have apparently constant thickness. 

Over the duration of the experiment, there is no evidence of any interfacial instability. 
These results are qualitatively similar to those of 42|,|43|], where C-940 solutions are displaced 
by glycerol solutions. The reason that there is no interfacial instability here is because the 
residual layers of C-940 solution are in fact fully static and unyielded. Such flows have 
been studied in some detail, theoretically and computationally 
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46 



471 ] . and are 



well understood. Even with a non-zero flow rate of C-940 solution, stable regimes may be 
predicted and found experimentally, see 48l. l49|. 

The two control sequences establish that without a chemical reaction and accompanying 
local rheology change, these displacement flows are stable. In other words, having a change 
in bulk rheology of the two fluids at an interface can result in instability but does not for 
the systems we study. 
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B. Chemically reactive unstable flows 

The flow behavior was significantly different from the control experiments in the reactive 
case, when the C-940 solution at pH = 3 was displaced by a saccharose solution at pH =11, 
at different flow rates. The initial interface penetrates in a sharp spike as before, but this is 
destabilized and the finger rapidly widens to nearly fill the pipe. A complex secondary flow 
develops at the interface between fluids. The flow seems to be dominated by large vortices 
advected by the flow, with a typical size of the order of the pipe radius. Typical images are 
shown in Fig. [6l 

As the front of the finger passes, the secondary flow instabilities persist along the sides 
of the finger. The secondary flow provides a feedback mechanism for the instability by 
bringing into contact new un-reacted fluid elements and taking away reacted highly viscous 
fluid. The initial pass of the finger front does not remove all the fluid 2 from the walls. 
However, the secondary flows result in a fairly rapid erosion of the residual layers. After 
the initial instability, small parcels of Fluid 2 pulled into the Fluid 1 stream react to form 
gelled solid regions that are advected along with the fluid. Close observation of video images 
reveals that some of these parcels appear to be in rigid motion, see e.g. Fig. & f . 

More detailed information on the structure of velocity field is obtained from DPIV 
images. We display in Fig. [7Ji & b two instantaneous velocity fields measured during a 
chemically reactive displacement experiment at roughly the midpoint of the fluid channel 
downstream. Before the entrance of the displacement front into the field of view, the flow 
is similar to Hagen-Poiseuille flow, see Fig. [TJa. The velocity profile is parabolic and there 
is no secondary fluid motion in a direction orthogonal to the mean flow. Analysis of time 
series of velocity fields (data not shown here) has shown that the time fluctuations are only 
due to instrumental noise, which accounts for less than 5% of the mean flow velocity. The 
structure of the flow field changes drastically after the passage of the displacement front 
through the field of view. As one can clearly see in Fig. [7b, the fluid motion follows a 
wavy-spiral pattern, with an apparent periodicity in the axial direction of roughly the tube 
diameter. The velocity field is now unstable and characterized by a rather strong radial 
component. Quasi-periodic entrance of slowly moving flow regions in the field of view can 
be associated with the passage of fully gelled solid parcels of fluid 2 illustrated in Fig. [6^ &f 
and discussed above. 
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Finally, we are able to confirm that the mixing process is locally effective and that the 
pH varies over some intermediate range in the displacement region. To do this we image 
unstable flows doped with pH sensitive colored dye, as shown below in Fig. [SJ Since the 
pif-dependence of the fluid rheology has been characterized, in principle we could use the 
DPIV images together with the data in Fig. [S] to construct local deviatoric stress fields. 
This technique however still requires some work. 

C. Displacement efficiency 

A more quantitative assessment of the displacement efficiency is obtained by processing 
LIF images to give an indication of the evolution in time of the finger width at a fixed 
position along channel, roughly 50 cm from the entrance. In each experiment a long time- 
series of images (~ 2000) is acquired, starting long before the entrance of the displacement 
front into the field of view and ending long after its passage. In these experiments fluid 2 
contained Polyamide spheres and fluid 1 a small amount of fluorescein (see section ITT1 C). 
The images prior to the entrance of the finger in the field of view are passed to the DPIV 
algorithm to obtain a time series of velocity fields. The time average of these images also 
provides a background image that is used to compensate for non uniform illumination of each 
image in the sequence. The fluorescent images of the displacement front/finger are digitally 
processed to extract the width (measured along the radial direction) of the finger. Each 
image is first converted to a binary image, and the edge of the interface is detected using the 
Sobel algorithm with a fixed brightness threshold. Spurious edge identifications, i.e. objects 
of the order of several pixels that are identified due to brightness inhomogeneities, are 
carefully removed using a morphological filter. Finally, the width of the finger is quantified 
from the width of the finger contour. 

During reactive experiments, the instantaneous shapes of the finger are usually non- 
symmetric in both azimuthal and radial directions. However statistically, (considered as 
an ensemble average over several positions of the finger), the symmetry does not seem to 
break. Thus, the finger width measurement can be interpreted as a volumetric measure of 
the displacement efficiency, although we prefer to display time series of the width. In Fig. [9] 
we display the time dependence of the width of the finger for several values of flow rate. 
For different flow rates the time until the displacement front enters the field of view will 
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be different. This residence time has been subtracted and the time (on the horizontal axis) 
is then scaled with the advection timescale, R/Uq. The finger width is normalized by the 
diameter of the pipe, d = 2R. 

Fig. [Hh shows the normalized width of the finger, Wj d, for the experiments in control 
sequence 1, where two Newtonian fluids are displaced. For each flow rate, W/d saturates at 
values smaller than unity, indicating that in all cases fluid 2 was only partially displaced. 
Qualitatively similar results are given in [ji]]. Fig. [9b shows analogous results for the reactive 
displacements at different flow rates. The interfacial instability clearly results in an efficient 
mass transport in the cross flow direction. For each flow rate, W/d saturates at values 
very close to unity. Observe that the initial points on the time series are noisy, followed 
by a rapid increase to near saturation, then slow approach to unity as the wall layers are 
consumed by the secondary flows/instability. The initial noisy part of the curves corresponds 
to destabilization of the initial interfacial spike. 

In comparing these two figures, recall that the time axis has been shifted to correspond 
to the initial appearance of fluid 1. In the case of the unstable displacements, this is harder 
to determine precisely due to the destabilized wispy spike. Once the bulk of the displacing 
finger arrives, the rapid increase to near the saturation values is similar in both cases. The 
wider spread of the curves in Fig. [9b is due to the early arrival of the initial spike and slower 
erosion of the final wall layer. In terms of a volumetric measure of efficiency, the reactive 
displacement leaves < 5% of the displaced fluid volume, and this appears to be insensitive 
to the flow rate. The stable Newtonian displacements leave between 20% and 60% of the 
displaced fluid in the pipe, with strong dependence on the flow rate. 

The key observations of Fig. [9b is that all the efficiency curves (time shifted for arrival) 
have similar shape and that the time scales with R/Uq. In terms of an averaged volumetric 
concentration of displaced fluid, say C(x, t) = (W/d) 2 , this suggests a functional form: 

C(x,t) = C(x - U t,tU /R), 

and the sigmoid shape is reminiscent of axial dispersion in a moving frame of reference. 
If this dispersion is governed by an axial diffusivity D a , the timescale for the spreading is 
t ~ R 2 1 ' D a . Thus, the collapse of our data with respect to the time variable R/Uq, suggests 
the scaling law 

Da ~ RU , 
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which infers a similar scaling law for a transverse diffusivity of mixing, D r , see §1 Al 

Even if redrawn volumetrically, the curves in Fig. [9b are not symmetric about the mid- 
concentration. Additionally, close inspection reveals small differences between the curves for 
different flow rates. These suggest that, although axial dispersion may provide a reasonable 
leading order model, the dispersion coefficient will vary nonlinearly with respect to C and 
have other weak dependencies. This is in fact physically obvious: the mechanism governing 
the initial instability of the wispy spike (C ~ 0) is clearly quite different to that governing 
erosion/break up of the residual wall layer of C-940 (C ~ 1). 

IV. HYDRODYNAMIC INSTABILITIES DRIVEN BY VISCOSITY GRADIENTS 

We would like to understand better the origin of the instabilities that develop into sec- 
ondary flows such as those in Fig. [6j Development of a complete stability analysis of the 
react ion- diffusion system coupled to the hydrodynamics is a formidable task, and perhaps 
unnecessary. One obstacle is that the reaction kinetics are not fully quantified. Another 
obstacle is to fully understand the coupling between the reaction and the hydrodynamics. 
This analysis is however underway and we hope to report the results later. For now, we offer 
some insights via analysis of a simpler toy problem. 

First, we note that the chemical reaction occurs on a time scale t c h which is much shorter 
than the hydrodynamic time scales of the problem. This is supported by the estimation 
below. 

Following ref. [sjj, the characteristic time scale at which the chemical reaction occurs 
is roughly t c h ~ 10 _9 s. In our experiments, the characteristic hydrodynamic times are 
the diffusion time, ip = R 2 /D, the advection time, t a = R/Uq, and the viscous time, 
t v = pR 2 /fio. If one considers D ~ 4.2 -10 -12 m 2 /s (which is roughly the diffusion coefficient of 
fluorescein in our solutions) and Uq ~ 5 mm/ s (which corresponds to half of our experimental 
range), one can estimate to ~ 7.2 • 10 6 s, t a ~ 1.1 s and t v ~ 0.3 s. The numerical estimates 
above show a clear separation of time scales in our problem: t c h *C t v < t a <C to- In view of 
this, it is reasonable to assume that the fast reaction instantaneously establishes an initial 
concentration profile when two fluids come into contact. This "reaction front" is a diffuse 
layer of reacted fluid separating the two bulk fluids, within which the pH is approximately 
neutral and the viscosity is consequently elevated. In the following we will investigate the 
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hydrodynamic stability of such a viscosity profile. 

This process is on the slow scale of t a . Note that the reaction is not auto-catalytic. Unless 
there is convective motion that brings fresh un-reacted fluids into contact, the reaction 
front broadens diffusively. No self-sustaining traveling chemical waves arise. In addition, 
unlike the many "chemical reaction plus buoyancy" studied referenced in ^ our fluids have 
matching densities that are unchanged by the reaction. Consequently, there is no momentum 
source resulting from the reaction. Therefore, the direct effect of the continual reaction is to 
create a source term in the advection-diffusion equation for the concentration. Indirectly this 
modifies the viscosity profile. For simplicity, however, we will ignore the effects of continual 
reaction, and focus on the fate of the initial, chemically created, concentration and viscosity 
profiles. 

Furthermore we ignore the pipe geometry 55j, and consider instead a symmetric plane 



Poiseuille flow in a channel of width 2R along which fluid is pumped at mean speed JJq. The 
fluid is assumed to have a concentration dependent viscosity 77(c) that consists of a base 
viscosity 170 that is augmented (by a factor a) over some finite range of concentrations, Aco, 
close to a fixed concentration value c = cq. For example, one such function would be: 

V(o)=vo[l + a cos [-^) , 

which is depicted in Fig. ITOk . Ignoring the reaction, the concentration satisfies an advection- 
diffusion equation and the fluid flow satisfies the Navier-Stokes equations. There exists a 
steady base flow in which the concentration varies linearly across the channel. The increased 
viscosity tends to flatten the base Poiseuille velocity profile within the diffuse layer, see 
e.g. Fig. dOb. 

We consider the linear stability of this base flow, using classical methods. The stability 
problem is governed by the Reynolds number Re = (pRUo)/f)Q < 1, the Schmidt number, 
Sc = fj /(pD) € [10 4 , 10 6 ], the initial concentration profile C , the layer thickness and posi- 
tion, Ac & c , and the amplitude a of the viscosity jump. Introducing the non-dimensional 
variables 

x y i U fj 

x = ^ V = "a> t = T , U = j-, v = — 
R R t v U Vo 

the growth or decay of a linear mode {f{y), c(y))e iax+at is governed by the following eigen- 
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value problem: 

a (D 2 -a 2 )f = ±-Lif + L 2 f + L 3 c (4) 
ac = iaRe {{DC Q )f - Ucj + — (D 2 - a 2 ) c (5) 

where D = d/dy, r\ = r)(C ), C = 1 — y, and 

Lx := -ia(u(D 2 - a 2 ) - (D 2 U)^j 

L 2 := (D 2 V ) (D 2 + a 2 ) + 2(Dr])D (D 2 - a 2 ) + rj (D 2 - a 2 ) 2 

We discretize the equations using Chebyschev polynomials and determine the maximal 
growth rates or as a function of wave number a (and the other dimensionless parameters), 



in the standard way, see e.g. 



sal- 



Examples of the results are shown below. 



A. Results 



Typically in our experiments we have Re < 1, so that we are far below the range for 
inertia-driven shear instabilities. If a = 0, the viscosity is constant so that equation (j3J) 
decouples and reduces to the classical Orr-Sommerfeld equation, which is stable at low Re. 
Equation (jSJ) then gives only stable modes with ctr < decreasing quadratically with a. As 
the amplitude a increases, unstable modes are found over an increasing range of a, including 
the long wave limit a = 0, with diffusion stabilizing the short wavelengths. Figure [TT] shows 
typical examples of the growth rate cjr vs a, for a range of different a, at two different values 
of Sc. 

The long wave limit a — > tends to give the largest growth rates, but these modes will 
not be excited in a developing finger-like displacement. Unstable wave numbers are found 
up to a critical wave number a c , which varies with a and Acq, increasing mildly with a. A 
possible physical interpretation of the (^-independent) long wavelength instability is simply 
that a different base flow may exist and be somehow more stable. In order to compare with 
the experimental results, note that the wave numbers are scaled with R~ l and the growth 
rates with the inverse viscous timescale, t v , (which is ~ 0.3 s for our experiments). Thus, we 
can see from Fig. [TT] that wavelengths ~ R are certainly excited at viscosity amplifications 
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a = 100. The growth rates are very rapid and we would therefore not expect to see linear 
modal effects. 

A complete exploration of the parameter space of this problem does not seem worthwhile 
and is not easy numerically. Although we typically have large Sc, the limit Sc — > oo results 
in loss of the diffusive terms in the concentration equation and it is these that stabilize the 
short wavelengths. The other limits of large a and small Acq will obviously eventually cause 
problems for the Chebyschev expansion, and would be better investigated analytically, as 
would the large and small wave number limits. 

V. DISCUSSION AND CONCLUSIONS 

In this study we have presented experimental evidence of an inertial free shear flow 
instability caused by changes in the local fluid rheology. We have focused on low Re number 
displacement flows in an horizontal pipe. An acid-base type chemical reaction occurs near 
the interface between fluids and results in molecular reorganization of the Carbopol-940 (C- 
940) polymer. The main result of this molecular reorganization is the formation of a stiff 
gel, characterized by a large viscosity and a significant yield stress. 

In order to test the effect of these rheological changes on the hydrodynamic stability of 
the flow, we have carried out two control experimental sequences, see Table U series 1 and 
2. In neither of the control experiments was a flow instability observed. If both fluids were 
Newtonian (first control sequence), the displacement flow was dominated by a long finger of 
fluid 1 penetrating into fluid 2, Fig. HI The flow was stable and the second fluid was never 
completely removed from the channel walls, Fig. [9|l. If the displaced fluid was a yield stress 
fluid, the flow remained stable at all times, (second control sequence), see Fig. |5] 

The behavior was significantly different in the reactive case. The interface between the 
fluids becomes unstable and mixes across the entire channel. The flow fields are unsteady 
and characterized by a strong secondary motion in a direction orthogonal to the mean flow 
direction, Fig. [7b. The relevance of our experimental findings is threefold. 

1. First, the coupling between chemical reaction and a strong non-monotonic local change 
in the fluid rheology is novel. For simple fluids, viscosity variations can be caused by 
temperature or concentration gradients, but these tend to be gradual and monotone. 
Use of complex fluids allows one to localize the change in viscosity and produce non- 
17 



monotone effects. 

The instability mechanism we have described in this paper is not restricted to interfa- 
cial flows, where acid-base type chemical reactions take place at the interface between 
two fluids via transport of unbalanced charges across the interface. Similar chemical 
reactions, which result in significant changes in the fluid rheology, may be triggered in 
a non invasive way, e.g. by exposing fluid parcels to either electromagnetic radiation 
or heat waves which will locally initiate a polymerization reaction. For example, there 
exist numerous polymeric materials (epoxies) widely used in photo lithography, whose 
viscosity increases dramatically (up to solidification, depending on the chemical nature 
of the polymer and the exposure conditions) upon exposure to monochromatic elec- 
tromagnetic radiation. Also, liquid silicone elastomers may display a similar behavior 
(they turn from viscous liquid state to solid elastic) when locally heated. This obser- 
vation, together with an instability mechanism similar to the one illustrated in this 
paper, may open new possibilities of externally controlling the hydrodynamic stability 
of inertia free shear flows. 

2. Secondly, we have demonstrated a significant increase in the displacement efficiency 
for reactive displacements. This has an impact on process applications such as these 
described in §1 A|. and this was the initial motivation for our study. 

3. Thirdly, we have demonstrated a new mechanism for low Re mixing that does not 
require large inertial energies. This in itself opens up many interesting avenues for 
practical applications, e.g. mixing in micro-fluidic flows. Indeed, efficient mixing in the 
absence of inertia requires an additional mechanism than molecular diffusion (which 
is the least efficient). Several micro mixing techniques have been proposed during 
the past half decade, all based on different mechanisms of generating secondary flows 
(steady or random in time). One of these techniques is based on the recently discovered 



elastic turbulence 



ex dynamic flow state 
"I 

541 ]. elastic turbulence 



52j . which is a random (in time) and comp 
in dilute solutions of linear flexible polymers. As shown in [53l . 
can be successfully employed in efficiently mixing viscous fluids in the absence of 
inertia. Although the flow configuration we used in this study is not a typical mixing 
configuration, based on the data presented in Figs. [6J and EJ one can suggest that 
the instability discussed in this paper could be alternatively used as a low Re mixing 



18 



method. 

Our hydro dynamic stability study is minimalist and we have presented only brief results. 
However, the simple toy model shows that long wavelength instabilities do result from non- 
monotonic viscosity variations of the magnitude that we have and in the correct range of 
Re and Sc. A more detailed study of these reactive displacements is underway. 
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FIG. 1: Schematic view of the experimental apparatus: FC-flow channel, GC-glass container, 
L-solid state laser, CO-cylindrical optics, LS-laser sheet, CCD -camera, Mi^-stepping motors, 
Vi 2,3-valves, TS-horizontal translational stage, SP-syringe pump. The lower inset illustrates 
schematically the flow configuration: the full curve represents the interface between fluids, the 
dotted lines indicate the reacting region and the horizontal arrows point the mean flow direction. 




FIG. 2: Shear rate dependence of the viscosities of the fluids: (o) displaced fluid 2 with pH ~ 3, 
(□) displacing Fluid 1 , (A) displaced Fluid 2 with pH 7. 
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FIG. 3: Dependence of the rheological properties of fluid 2 on pH: (a) viscosity measured at 
7 = 2s _1 (b) yield stress. 






FIG. 4: Example fluorescent images of the interface in an experiment from control sequence 1: (a- 
b) Fluid 1 - 65% saccharose solution, Fluid 2 - 66% saccharose solution, flow rate Q = 0.145 ml/s. 
The two images are separated in time by 5 s. 
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FIG. 5: Example images of the interface in an experiment from control sequence 2 (Q = 0.08 ml/s). 
Fluid 1 has been doped with fluorescein for visualization of the fluids interface. 
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FIG. 6: (a-d) Fluorescent images of the interface in a reactive displacement: Fluid 1 - 65% saccha- 
rose solution, Fluid 2 - 0.1% C-940 in 66% saccharose solution, (e-f) Fluorescent flow images long 
after the entrance of the unstable interface in the field of view; the images are separated in time 
by approximately 5 s. The dotted lines highlight gelled structures tumbling downstream. 
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FIG. 7: Instantaneous velocity fields for a chemically reactive displacement experiment, Re ~ 
0.15: (a) before entrance of the finger in the field of view (b) after the entrance of the finger in the 
field of view. The false color maps represent the magnitude of the velocity. In order to enhance 
the clarity, we only plot half of the vectors and re-scale the length of the arrows. 




FIG. 8: Iso pH regions near the interface: yellow (region l)-pH ~ 3, green (region 2)-pH ~ 7, 
blue (region 3)-pH « 11. 
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FIG. 9: (a) Normalized width of the tip versus the normalized displacement distance, Uot/R, for 
several values of the flow rate: (■) Q = 0.063 ml/s, (•) Q = 0.145 ml/s, (A) Q = 0.19 ml/s, (<) 
Q = 0.3 ml/s. The experiments belong to Control sequence 1 (see table [TJ) . (b) Normalized width 
of the tip versus the normalized displacement distance, Uot/R, for several values of the flow rate: 
(<) Q = 0.13 ml/s, (■) Q = 0.18 ml/s, (▼) Q = 0.2 ml/s, (•) Q = 0.31 ml/s, (►) Q = 0.47 ml/s. 
The experiments belong to the Reactive sequence (see Table H]). 
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FIG. 10: Example basic viscosity and velocity profiles for cq = 0.5, Acq = 0.1 and a = 100: (a) 
Normalized base viscosity profile, (b) Normalized base velocity profile. 
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FIG. 11: Example results: ctr vs a for the base flow of Fig. [10] with the viscosity amplitudes 
(A) a = 100 (■) a = 5 and (•) a = 20: (a) Sc = 10 4 . In the inset we display the dependence of 
the critical wave number, a c , on the amplitude a; (b) Sc = 10 6 . 



Tables 





Fluid 1 






Fluid 2 










Composition: 




Composition: 










65 % saccharose 




66 % saccharose 








Exp. 


r)\ = 114 mPas 














Sequence 


pH 


Re 


C-940 (%) 


pH 


Re 


m (Pas) 


Q (ml/s) 


Stab. 


Control 1 


7 


0.03 - 0.18 





7 


0.03-0.15 


0.138 


0.06 - 0.3 


S 


Control 2 


7 


0.02-0.3 


0.1 


7 


(0.1 - 1.5) • 10" 03 


20 


0.03 - 0.4 


s 


Reactive 


7.7-11.5 


0.07-0.2 


0.1 


3 


0.05-0.3 


0.144 


0.1 - 0.5 


u 



TABLE I: The experiments and the corresponding fluids properties. For each fluid, Re was calcu- 
lated using the DPIV measured mean flow velocity, the corresponding viscosity coefficients and 
the densities p\ ~ p2 ~ 1285 kg/m 3 . The flow rate Q has been estimated via the DPIV measured 
mean flow velocity. In the last column of the table, the symbols S, U stand for stable and unstable, 
respectively. 
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